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Abstract 

We study a d-dimensional lattice model of diffusing coalescing massive particles, with 
two parameters controlling deposition and evaporation of monomers. The unique stationary 
distribution for the system exhibits a phase transition in all dimensions d > 1 between a 
growing phase, in which the expected mass is infinite at each site, and an exponential phase 
in which the expected mass is finite. We establish rigorous upper and lower bounds on the 
critical curve describing the phase transition for this system, and some asymptotics for large 
or small deposition rates. 
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1 Introduction 

In section 1.1 we state the rules for the interacting particle system studied throughout the 
paper, describe the main new results on the phase transition, and explain the intuition behind 
the results. In section 1.2 we describe some previous work on this and related models, including 
numerical and theoretical intuition into the two phases. 

1.1 Summary of main results 

Particle Rules. Particles live on Zi d . Each particle has a mass with value in N = {1,2,...}. 
Particles move between nearest neighbour sites as independent rate one simple random walks. 
When a particle moves onto an occupied site, the two particles instantly coalesce producing 
a single particle whose mass is the sum of the masses of the two interacting particles. Thus 
each site contains either zero or one particle at all times. There are two additional dynamics 
controlled by parameters p, q > 0. 

Evaporation. Rate p evaporation, independently for each particle, reduces their mass by 

one. A particle of mass one therefore disappears at the time of an evaporation. 

Deposition. At rate q, independently at each site x G Z d , particles of mass one, often called 

monomers, are deposited. A monomer deposited onto an occupied site will instantly coalesce, 
increasing the mass of the existing particle by one. 
There is a unique stationary distribution for this process, and the law of any solution converges 

to this stationary distribution (see Proposition 3). However the properties of the stationary 



distribution depend strongly on the values of the parameters (p, q) . In the next section we 
summarize some earlier work which describes two phases for this stationary distribution, called 
the exponential phase and the growing phase. In this paper we characterize these phases by the 
expected mass at a site x € f L d . Note the stationary distribution is translation invariant. 

Theorem 1. Let M^O) have the distribution of the mass of the particle at the origin in the 
stationary distribution. Then, there is a non- decreasing function q c : [0, oo) — > (0, oo) so that if 
q G [0,q c (p)) then E[Moo{0)] < oo, while if q > q c (p) then E[M oo(0)] = oo. 

The existence and non-triviality of q c (p) will be shown, in section 2, by analysing the equa- 
tions for the first and second moments of the mass at a fixed site. The moment equations are 
not closed. However by exploiting global properties of the distribution, namely monotonicity 
and pairwise negative correlation, they yield a differential inequality for the first moment that 
implies the non-triviality of q c {p). This was a surprise to the authors, and it seems to be the 
presence of negative correlations, natural for systems undergoing coalescence, that allows the 
moment equations to be replaced by differential inequalities that are in the correct direction 
to yield non-trivial information about the phase curve. Contrast this with the contact process, 
which has positive correlations and where moments equations do not seem to help in establishing 
the phase transition. For large or small values of the parameters p, q the system should simplify. 
In section 3 we give quite detailed theoretical reasoning behind the following asymptotics. 

Conjecture. Large p,q asympotics. In all dimensions q c (p) < P and 

Urn P ~^ P) =&(d)€(0,oo). (1) 

Small p, q asympotics. In dimension d = 1 

lim q c (p)p~ 3 ^ 2 = etc S (0, oo), (2) 
p^O 

and in dimensions d > 3 

\\mq c (p)p- 2 = {Ap d )- 1 (3) 
p^O 

where pd is the escape probability for simple random walk on Z d , that is the probability that 
a walk leaving the origin never returns to the origin. For each of the three limits, there is 
a continuum approximation that suggests the answer. The intuition behind the asymptotic 
for large p, q is that the many mass changes occur between each random walk step, so that the 
model can be approximated, on suitable time-mass scales, by a simpler model with continuously 
evolving masses. We will show that this simpler model has a one-parameter phase transition 
with critical value /3 c (d). For small p, q in d = 1 there are many random walk steps between 
each deposition and the model can be approximated, on suitable space-time scales, by one with 
coalescing massive Brownian particles, and this model too has a one-parameter phase transition 
with critical value a c . For small p, q in transient dimensions d > 3, where there is no continuum 
approximation for the dynamics. However the when p is small and q = ap 2 the occupation 
density in the stationary distribution is O(p). By rescaling the stationary distribution in space 
by p l / d we expect a continuum approximation which has a compound Poisson distribution, 
that is particles are positioned at a Poisson rate s a , and the particles have independent masses 
attached with a law determined by a generating function <p a . There is a phase transition for 
these compound Poisson limits which can be analysed, following the route used in [12] for the 
mean field analysis of this model. The critical value for the continuum approximation can be 
exactly found as a = (4j>d) _1 , suggesting the asymptotic (3). The idea that there are various 
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one parameter models that reflect the limiting situations in the two parameter phase diagram 
is similar to the situation for a branching reacting system studied in [14]. In this paper we 
do not establish all the weak convergence arguments needed to establish the conjectured exact 
asymptotics, but the moment methods are already sufficient to imply the following weaker forms 
of the asymptotics for some cases: 

Theorem 2. 

< hmmt — < hmsup -j- — < oo, m all a, (4) 

p^oo p 1 /^ p^Qo pW^ 

< liminf q c (p)p 3 ^ 2 < Hmsup q c (p)p~ 3 ^ 2 < oo, ind=l. (5) 
1.2 Background for model 

An earlier paper on this model [6] gives a brief bibliography for models of various physical 
phenomena that involve diffusion, coalescence (or aggregation) and deposition (or immigration) 
including the survey [13]. The river network model in [19] and the Takayasu model [20] are close 
to our model but with zero evaporation. This zero evaporation model is well understood and has 
a (non-reversible) stationary distribution with a constant mass flux from small to large masses. 
The mass distribution at a fixed site has polynomial tail, which behaves like m -4 / 3 for large m 
in d = 1 and m _3//2 in d > 3. See [17, 4] for recent developments. The effect of evaporation 
was studied in [12], which demonstrated numerically the phase transition that we study in this 
paper and included a mean field analysis. A more detailed description of the two phases is 
given in [6], [5] where the mass balance in stationarity is analysed formally, and confirmed by 
numerical investigations. In the exponential phase, where q < q c (p), the mass distribution at a 
site has exponential tails. In the growing phase, where q > q c {p) the mass distribution at a site 
has polynomial tail and mimics the tails when there is no evaporation. Indeed it is conjectured 
that the entire space-mass distribution, at large masses, should be well approximated by that 
of the model with zero evaporation but with a modified deposition rate q' = q — ps(oo), where 
s(oo) is the occupation probability for any site in the stationary distribution. Note that in [6] 
the phase transition is considered as q c (p) = mf{q : s(oo) < q/p}. (One can show, see appendix 
4.1, that s(oo) < q/p for s > q c {p) so that this defines a true transition). If s(oo) < q/p then 
the first moment equation shows that the first moment is infinite in the stationary distribution, 
and hence q c (p) < <7c(p)- We do not know whether q c (p) = q c (p), and we believe this might be 
a useful tool in proving some of the many predictions about this model: exponential moments 
throughout the exponential phase, behaviour at criticality (see the scaling predictions in [6]), 
regularity of the critical curve p — ^ Qc(p) e.t.c. Numerical simulations shown in Fig. 1 illustrate 
the relation between our bounds and the true phase boundary in d = 1, 2, 3. It is often observed 
that mean field approximations become exact in the limit d — > oo and so one may conjecture that 
the critical curve approaches the lower bound (13) as d grows. An approximate d = oo model 
has been simulated using particles that jump to any other lattice site with equal probability, 
and the simulated critical curve for this model supports the large d conjecture. 

2 The phase diagram 

In section 2.1 we state monotonicity properties, which follow from standard methods as in [10], 
and negative correlation property for the model. In section 2.2 we use these and the moment 
equations to show the existence of the critical curve q c (p)- The proof of negative correlation, 
based on the BKR inequality as in [2], is given in section 2.3. 
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Figure 1: Numerical investigation of the phase boundary. Squares, white circles and black circles 
are numerical measurements of q c (p) in d = 1,2,3 correspondingly. Triangles are numerical 
measurements of phase boundary in the approximate 'd = oo' model. The solid curve is the 
rigorous lower bound. 

2.1 Tools 

We write Mt(x) for the mass of the particle at site x at time t. The construction of the process 
Mt = (Mt(x) : x G Z rf ) started from suitable initial conditions follows from standard arguments 
on particle models (see for example chapter IX of [10]). One method is to consider M as the 
solution of the following system of stochastic differential equations driven by Poisson processes 
and indexed over the lattice: almost surely, for all t > 0, x G Z d 

dM t (x)=dP t {q \x)-x(M t _(x) > d)dP { t p \x)-Y J M t -{x)dP t (x,y) + Y J M^(y)dP t {y,x) (6) 

y~x y~x 

where y ~ x means that y and x are nearest neighbours, and the independent families of Poisson 
processes (P^ q \x) : x G Z d ) of rate q control the deposition, (P^^x) : x G Z d ) of rate p control 
evaporation and (P(x, y) : x, y G Z d ) of rate (2<i) _1 control the random walk steps. We consider 
only solutions where for all x the path t — > Mt(x) is cadlag (right continuous with left limits). 
Furthermore solutions should be adapted to a filtration (T t ) where, for any s < t, the increment 
Pt — P s of any driving Poisson processes is independent of J 7 ,. Such solutions can be found 
without any further restrictions on the initial conditions. Indeed we may take as state space 
the product space S = N z (with the product topology). The point is that for independent 
random walkers without coalescence there would normally be a growth condition at infinity on 
the initial state, but this is unnecessary for instantly coalescing particles (due to the bounded 
jump intensities). We collect together some basic results on existence, uniqueness, stationary 
distributions, moments, monotonicity and correlation, giving comments on the proofs at the 
end of this section. 
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Proposition 3. (i) For any initial condition Mq G S that is independent of the driving 
Poisson processes, there exists a pathwise unique solution to (6) with initial condition Mq. 
The laws of solutions form a Markov family with a Feller semigroup. 

(ii) If for some 6 > and k G N the initial moments E J2 X \Mo(x)\ k exp(— 0\x\) are finite, 
then one has finite moments at all times 

£supV \M t (x)\ k e~ 9 ^ < oo, for all T. (7) 

(Hi) There exists a unique stationary measure on S. For any solution as above, the law of M t 
converges in distribution, as t — > oo, to the stationary distribution. 

A stochastic monotonicity property holds, frequently used for interacting particle systems, 
as follows. Put a partial order on S by writing rj < fj for if rj(x) < fj(x) for all x G Z d . A 
function F : S — > R is called non-decreasing if F(n) < Fiji) whenever rj <rj. 

Lemma 4. Let (M t : t > 0) be the solution to (6) started from Mq = 0. Then for all measurable 
and non- decreasing F : S — > [0, oo) the expectation E [F(M t )] is non- decreasing in t and q and 
non-increasing in p. 

Lemma 5. Let (M t : t > 0) be the solution to (6) started from a deterministic initial condition 
Mq = rj G S. Then for any t > the variables (Mt(x) : x G Z d ) are negatively associated. In 
particular, 

E \f(M t (x))g(M t {y))) < E [f(M t (x))] E [g(M t (y))} (8) 
for measurable and increasing /, g : N — >■ [0, oo). 

Comments. 1. Proposition 3 and the monotonicity Lemma 4 are rather standard results for 
(suitable) particle systems (for example see chapter 2 of Liggett [10]). We give some details in 
appendix 4.1, where we emphasize the use of differential equation comparison methods. Indeed 
there is a basic comparison theorem for solutions to (6): if two solutions (Mt : t > 0) and (Mt : 
t > 0) satisfy Mq < Mq a.s., then Mt < Mt for all t > a.s. This can be obtained by standard 
differential equation methods, namely a Gronwall estimate on E^2 X e~ e ^x(Mt(x) > Mt(x)). 
A similar comparison theorem, with suitably coupled Poisson drivers, yields the monotonicity 
in p, q in Lemma 4. The monotonicity in t comes from the fact that Mq = is a minimal initial 
condition. 

2. The monotonicity can be used to show the existence of a stationary distribution, again in 
familiar way for attractive particle systems. Convergence starting from zero initial condition 
follows from the monotonicity of the Laplace transforms £'[exp(— J2 X Mt(x)4>(x))], for 4> : Z d — > 
[0, oo) of compact support. One needs however to check that the limit variables (M OQ (x) : x G 
Z d ) are non-degenerate. Informally, if P[M 00 (x) = oo] > then there would be a steady state 
for an extended particle system that allowed particles with infinite mass. But in this extended 
system the infinite mass particles do not feel collisions with finite mass particles, and hence 
act as an autonomous coalescing system, with no immigration. Such a system does not have 
a non-zero steady state. Pushing this argument further, one can find a maximal entrance law 
for the extended system where all sites start with infinite mass. The decay of the infinite mass 
particles leads to the coupling of the maximal and minimal entrance laws and this implies the 
facts about the stationary distribution. Details are in the appendix. 

3. Negative association is the strongest of various negative correlation type properties for 
random vectors, with pairwise negative correlation (8) being perhaps the weakest. Pairwise 
negative correlation, which is all we need for our arguments, was established for coalescing 
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random walks in Arratia [1] by a Markov duality argument. We will prove Lemma 5 by a small 
modification of the arguments in van den Berg and Kesten [2] , which deals with purely coalescing 
system, and which uses the van den Berg, Kesten, Reimers (BKR) inequality. Although their 
argument goes through without much trouble for our model, we do have to check that the 
extra deposition and evaporation mechanisms do not destroy the proof. This is a key tool for 
this paper, and rather less standard, and so we give the proof in section 2.3. Some restriction 
on initial conditions is necessary. The examples in Liggett [11] show that negative correlation 
properties do not mesh well with random initial conditions, which is one reason why, unlike 
positive correlations, differential techniques have not been that successful in establishing such 
properties. 



2.2 Phase transition via moments equations 

Throughout this section we consider zero initial conditions. The distribution of (M t (x) : x € Z d ) 
at each fixed t is then translation invariant on Z d (by the uniqueness of solutions to (6)). 
Monotonicity, as in Lemma 4, implies that the moments E[M[ r(0)], for k > 0, are non-decreasing 
in t and q, and non-increasing in p. Set 



q c (p) =inf \q : lim E[M t {0)] 



OO 



where the expectation uses zero initial conditions and parameter values (p, q). This defines q c (p) 
as a non-decreasing function with values in [0, oo]. Once we have shown that q c (p) £ (0, oo) the 
remaining statement in Theorem 1 follows from monotonicity. We use the differential equations 
for the first moment m\(t) = E[M t (0)] and the second moment rri2{t) = E[M^(0)]. These can be 
derived from (6) by developing M t fc (0) using calculus and then taking expectation and exploiting 
translation invariance. They imply that the first moment m\{t) is continuously differentiable 
and satisfies 

-±(t)=q-ps{t) (9) 

where s(t) = P[M t (0) > 0], the probability the origin is occupied at time t. In particular 
if q > p then mi(t) f oo showing that the steady state will be in the growing phase where 
i?[Moo(0)] = oo. This shows q c (p) < p, but it is simple to improve this upper bound. Indeed 

ds 

-(t) = q{l - s(t)) - pP[M t (0) = 1] + P[M t (0) = 0, M t (e) > 0] - s(t) 

where the last two terms arise from the simple random walking into and out of the origin, and 
e denotes a nearest neighbour to the origin. Bounding 

P[M t (0) = 0, Mt(e) > 0] < P[M t (0) = 0] = 1 - s(t) 

and discarding the term involving p one reaches 

-(t)<(l+q)-(2 + q)s(t). 

Since s(0) = this implies that s(t) < (1 + q)/(2 + q) for all t > and substituting this into 
the first moment equation (9) one sees that mi(t) t oo whenever q(2 + q) > p(l + q). This in 
turn shows that 

( P -2) + v^r4 

9c(p) < o • (10) 
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Further incremental improvements like this are possible on the upper bound, but without any 
particular hope that one is approaching the true critical curve. For a lower bound we start with 
the second moment equation 

^(t) = g(2mi(t) + 1) -p(2mi(t) - s(t)) + 2E[M t (0)M t (e)\ (11) 

where the last term arises from the particles leaving or entering the origin. By monotonicity 
dm,2/dt(t) > and by negative correlation E[M t (0) M t (e)} < m\{t). With these substitutions 
we find 

2m? (t) - 2(p- g)mi(t) +q + ps{t) > 0. 

The first moment equation (9) and monotonicity of mi(t) imply that s(t) < q/p for all t. This 
last substitution yields a quadratic inequality for the first moment: 

m\{t) - (p-g)mi(t) + q > for alH > 0. (12) 

Note that the initial condition is mi(0) = and that t — )■ mi(i) is continuous and non- 
decreasing. Therefore whenever the quadratic in (12) attains a strictly negative minimum 
on (0, oo) the first moment mi stays bounded by the first positive real root for all t > 0. This 
happens precisely if p > q and (p — q) 2 > 4q, and the corresponding upper bound is 

mi (t)< {p - q) -^ { P- q)2 - Aq foralU>0. 

Resolving the condition (p — q) 2 > 4q with respect to q < p we find that when q < p+2 — 2y/p~+l 
the first moment stays bounded and hence 

9c(p) > P + 2 - 2v^TT. (13) 

Note that this analysis is dimension independent. Moreover the formula for the lower bound 
appears in the mean field analysis performed in [12]. This is no surprise since the negative 
correlation allows us to replace the second moment equation by a differential inequality that 
agrees with the mean field approximation. It is natural to try to apply the same differential 
inequality methods to investigate other moments, but we have not profited much from this. 
Fro example, exponential moments E[e 9Mac ^] are conjectured to be finite, for suitable 9 > 
0, whenever q < q c {p)- Consider the following exponential moment, under the zero initial 
condition, for a fixed z > 1. 

i/>(z,t) = E[z Mt ^ - 1] 



This solves formally (since the finiteness of such moments is not yet established) the equation 

^(t) = -(pz^-q^z-^m + il-ps^z-^z-l) 
+ E[z Mt{0)+Mt{e) - z Mt{0) - z Mt(e) + 1]. 

Monotonicity implies that ip(t) is increasing and negative correlation shows that 

E[z M t (0)+M t (e) _ z M t (0) _ z M t (e) + ^ < 

Using this we obtain the following inequality involving the exponential moment ip(t) and the 
non-occupation probability s(t): 

ip{t) 2 - (pz" 1 -q)(z- l)ip(t) + (q- ps{t)z~ l ){z - 1) > for all t > 0. 
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Now the previous upper bound on s(t) is acting in the wrong direction, but using s(t) > we 
reach the quadratic ip(t) 2 — (pz^ 1 — q)(z — l)tp(t) + q(z — 1) > 0. As before, if this quadratic in 
ip has a strictly negative minimum at a positive value of ip then the exponential moment stays 
bounded. It is easy to see that this holds when q > is small, for suitable z G (l,q/p). This 
suggests that the method could be used to show some information on exponential moments, 
but only for q far from the critical curve. 



2.3 Proof of negative association 

Various notions of negative correlation are studied, and they have been exploited for a variety of 
interacting particle systems. See Newman [15], Pemantle [16], Liggett [11] for overviews. Real 
variables (X±, . . . , X n ) are called negatively associated if for any two disjoint subsets Ji,J 2 ^ 
{1, . . . , n} and any non-decreasing fi : R Ji — > [0, 00) and one has 

E[h{X 3 -j G J 1 )f 2 (X j : j G J 2 )] < ElMXfJ G Jx)]E[f 2 (Xj : j G J 2 )]. 

Infinite vectors of variables are called negatively associated if each finite subset is negatively 
associated. We give here the details for the proof of Lemma 5, following the strategy from 
Lemmas 2.4-2.7 in van den Berg and Kesten [2]. They exploit the BKR inequality, an inequality 
for product measures on finite sets, which we state here. Suppose V, S are finite sets and that 
/i is a product measure on the product space S v . For K C V and ui = (uj v : v G V) G S v 
define [lo]k = G S v : uj v = oj v for v G K}. The set [oj]k is called a cylinder set with base set 
K. For A, B C Q define ADB to be those u for which there exist disjoint K , K' C V so that 
[uj] k C A and [cu] K > C B. The BKR inequality states (see [18]) that n[A\JB] < fx[A]n[B]. Fix 
t > 0. A graphical construction for the process over [0, t] can be given using the same Poisson 
processes that drive the equation (6). This graphical construction contains, in a direct fashion, 
information on the genealogy of individual particles, and we recap the required notation. For 
each x ~ y, at each jump of P(x, y), place on the space time lattice [0, t] x Z d an arrow starting 
at (t,x) and ending at (t,y). This arrow corresponds to the particle at (t—,x), if it exists, 
jumping to site (t,y). Using only a realization of these arrows there is a natural notion of a 
path from (s, x) — > (s',x') when < s < s' < t. If such a path exists a particle at (s,x) 
would end up at (s',x') (ignoring the mass labels for the moment). See, for example, chapter 
3 of Durrett [7] for a careful definition. Using the other Poisson drivers (p( g \x) : x G Z d ), 
(p(p)(a;) : x G Z d ) we add Poisson points where evaporations and depositions occur. Almost 
surely, one can then trace through the evolution of the masses (rj(x) : x G Z d ) at time zero 
to yield the final masses (Mt(x) : x G Z d ). We define an embedded discrete time structure 
which will yield the finite product structure required for BKR. Fix N,L G N. Define variables 
(N(k,x), Z(k,x) : < k < N - l,x G Z d ) as follows: let J fcjJV be the time interval (f , ^^]; 
write P((r, s]) for the increment P s — P r ; let 

N(k,x) = ^P(x,y)(I k ,n) + P^ P \ x )(h,N) + P ( ' 9 Hx)(Ik,N), 
namely the number of Poisson events that occur at x during the time interval 1^ jy; and let 



Z(k, x) 



if N(k,x) = 0, 

p HN(k,x) = landP ( P\x)(I k>N ) = l, 

q if N(k,x) = 1 and P(i\x)(I k ] N ) = 1, 

±ej if N(k, x) = 1 and P(x, x ± ei)(I kj N) = 1 ; 

A if JV(Jfe,ar) > 2, 
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where (±ej : i = l,...,d) are the unit vectors in Z d . Let S be the finite set of labels 
{0,p,q,±ei, A}. Let V = {0,...,N- 1} x (Z d n [-L, L] d ). Then the vector 

Z := (Z(k, x) : < k < N - 1, \x\ < L) 

has the desired product law on S v . Fix finite disjoint sets J\, J2 C Z d . Choose Lo so that 
\x\ < Lq for all x G Ji U J2. Fix non-decreasing fa : R J ^ — > [0, 00), and cij > 0, for i = 1, 2, and 
set 

4 = {/i(M t (xj) : j G Ji) > a,}, for i = 1, 2. 

Negative association follows if we can show P[A\ n ^b] < -P^iJ-Pl^b]- We define a good set 
Q = Q l L N n Q\ as follows. 

^i,AT = \ N(k, x) + Y, N (k, y) G {0, 1} for all (k, x) G V 1 

is the set where at most one Poisson event occurs during each interval Ik,N for each \x\ < L, and 
furthermore if one occurs then neighbouring sites have zero Poisson events during this interval. 

Q L = {there is no path (x,s) — > (y,t) for any s G [0,i), \y\ < Lq and |x| > L} 

is the set where no particle moves from outside [-L, L] d to affect the values of (M t (y) : \y\ < Lq). 

Claim 1. (i) P[Q 2 L ] -> 1 as L -»• 00; (ii) P[<?i )JV ] -»• 1 as N -»• 00. Write (Xf : t > 0) for 
a simple rate one random walk started at x. To estimate it is enough, for each x with 

\x\ > L, to consider only paths started at the points (x, 0) or (x, s) where s < t is a jump time 
of P(x, y). The expected number of such points that start paths leading to a (y, t) with \y\ < Lq 
is bounded by 

£(l + t)supP[Xj€[-Lo,Lo] d ] 

|x|>L ^' 

and simple random walk estimates show this approaches zero as L grows. Part (ii) is a simple 
Poisson calculation. Let Qq be the (good) subset of S v where the value A is never taken 

and where no pair of neighbouring sites, (k,x), (k,y) with x ~ y, simultaneously take values 
different from 0. On Q\ N we know that Z G fio- 

Claim 2. There exist A 1: A 2 C O C S" 17 so that 

AidQ = {Z G ^4j} fl fori = 1,2. 

Indeed, we may take = {oj G Slo : = oj} n ^ C n ^}. To see this, it is sufficient check 
that on the set Q the values of (Mt(x) : |x| < Lq) can be reconstructed only from Z and the 
the initial masses (r/(x) : \x\ < L). Of course, this is the purpose of the good set. For oj G 0,q 
we may define (inductively in k) mass values (m(k, x,oj) : \x\ < L, k = 0, . . . , N — 1), starting 
at k = with the initial masses (ry(x) : |x| < L). This is the natural discrete analogue of our 
process. If oj^^x takes the value p or q the value m(k — 1, x, w) can be updated to m(fc, x, w) in the 
appropriate manner (since no neighbouring site will affect the mass at x during this interval). 
If ujk^x = iej the mass m(k,x,uj) is set to 0, and the mass m(& — l,x,w) is added to the mass 
at y = x ± ej (if |y| < L). Again the restriction to ensures there is never a conflict due to a 
pair of events occuring which would require us to decide which action (deposition, evaporation, 
movement) to do first. The mass values m(k,x,uj) do not take into account any mass entering 
from sites \y\ > L. But on the set Q\ we may ignore all Poisson points at sites x G" [—L,L] 2 



9 



without affecting the values of (M t (x) : \x\ < Lq). Hence the final values (m(N, x, Z) : \x\ < Lq) 
will agree, on Q, with (Mt(x) : \x\ < Lq). 

Claim 3. There are maps K- L : Slo — > V(V) (where V(V) is the power set of V) so that 

M= |J [u\ Ki{w) fori = 1,2, (14) 

and 

K 1 (u)nK 2 {u) =0 for all u G S v . (15) 

There is a natural notion of paths on £Iq. For a fixed w £ Oo, and for < k < k' < N and 
|x|, \x'\ < L, we write (k, x) (k', x') if there exist (zj : j = 0, . . . , n) G Z d n [— L, L] d satisfying 

(1) zo = x, z n = x' , n = k' — k, 

(2) \xj -Xj-!\ G {0,1} for j = l,...,n, 

(3) if - Xj_i = dzej then uj k+j ^ Xj = ±e, for all j = 1, . . . , n, 

(4) if - Xj-i = then uj k+ j, Xj G {0,p, <?} for all j = 1, . . . , n. 

In words, the values of Lo k+ j iZj along j = 0, . . . , n guarantee that a particle starting at (k, x) 
would end up at (k',x'). Sites (k,x) can have a (unique) path to at most one site (N,y) (some 
will exit the region [— L,L] d ). Any two such paths will coalesce at the first site in common. 
Now define, for uj G 0.$, 

Ki(ui) = {(k, x) : (k, x) ~~> (TV, y) for some y G Jj} . 

Then Ki(oj) is the union of the coalescing paths that lead to a site in Jj. Claim (15) is immediate 
since a path cannot lead both to J\ and J 2 . It remains to check that if uj G Ai then C A{. 

Observe that (m(k, x,oj) : x G Jj) can be calculated using only on the values (uj ktX '■ (k, x) G 
Ki(u})), indeed m(N, y, uj) is found by tracking the mass of the initial condition and depositions 
along paths that lead to (N,y). For uj' G we have Ki(uj) C Ki{uj') by definition of a 

cylinder set, with a strict inclusion possible if other sites are now connected to ((TV, y) : y G Jj). 
But these extra sites can only lead to larger values of mass ending up in Jj, that is m(N, y, uj') > 
m(N,y,uj) for all y G J\. Since fi are non-decreasing we must have 

fi (m(N, y, J) :y G J,) > /j (m(N, y,u) : y G J,) > a, 

and hence {2 = w'} n 6 <^ H <5 as desired. 

Proof of negative association. We proceed as in Lemma 2.4 of [2]. Using claims 2 and 3 

A 1 n A 2 n g c 2eAini 2 

= z G J MtfiM n (J M/^M 

cjGAi UJ&A2 

C Z G AiDA 2 , 

The last line following from (15). Hence, applying BKR to the law of Z, 

P[iiDl 2 ] < G AiDA 2 ] + p[g c ] 

< P[i? G A\]P[Z G -A 2 ] + P[£ c ] 

< P[i!]P[i 2 ]+3P[S c ]. 

Now claim 1 completes the proof. 
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Remark. One might start by proving negative association for a purely discrete time process 
on a finite lattice, denned on a finite product structure, and this avoids the need for the good 
set Q. One then might show that the such discrete time processes converge to the desired 
continuous limit, noting the conclusion carries over naturally under convergence in distribution. 
The trick of embedding a discrete time structure inside a graphical structure for the desired 
infinite lattice continuous process, and the good set Q, is designed specifically to avoid a lengthy 
weak convergence argument. 

3 Asymptotics 

Examining the upper bound (10) and lower bound (13) on q c (p) for large and small p, q one 
finds the asymptotic inequalities (that is in terms of limiting quotients) 

V - 2y/p ^ q c (p) r< P - 1 as p ->■ oo, 
i < q c (p) <l asp^O. 

In this section we show that more accurate, dimension dependent, asymptotic information on 
the critical curve is often possible. 

3.1 Large p, q 

For large q,p there are many evaporations and depositions between each random walk step. 
The correct approximating model is one with continuous masses in [0, oo) at discrete lattice 
sites x G Z d . We first study this continuous mass model. 

3.1.1 Continuous mass model. 

Masses (X t (x) : x G Z d ,t > 0), indexed over sites in Z d , evolve as follows. Between jumps the 
masses X t (x) evolve according to the stochastic differential equations 

dX(x) = -(3dt + V2dB(x) + dL(x), (16) 

driven by independent Brownian motions (B t (x) : x G Z d ,t > 0). The parameter (3 will reflect 
the excess of the evaporation rate q over the deposition rate p in the original model. The term 
L{x) is the local time of X(x) at zero, so that dL{x) charges only the set {t : X t (x) = 0} and 
ensures that Xt(x) > 0. Additionally the mass at each site x has independent rate one Poisson 
driven jumps, where at a jump time t a nearest neighbour y is chosen at random and the masses 
change to 

X t (x) = 0, X t (y) = X t _(x) + X t _(y), 

that is the mass at x jumps onto y and coalesces. We now argue that the moment method 
again establishes the existence of a phase transition with a critical parameter (3 c {d). Since 
we use it only as a guide to a later calculation for the discrete mass model, we do not give a 
careful construction of the model or prove the required monotonicity and dependence properties 
(although one way to establish them is to take limits in the corresponding properties for discrete 
mass model). Define moments for the system, with zero initial condition, by rrik(t) = E[X^(0)] 
and set the critical value as 

p c {d) = inf{/3 > : mi(oo) = oo}. 
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We aim to show f3 c (d) G (0, oo). The second moment equation works smoothly. Ito calculus 
and negative correlation show that 

= 2 - 2pm 1 {t) + 2E[X t (0)X t (l)} < 2 - 20mi + 2m 2 . 

The argument from section 2.2 shows that m\ stays bounded if (3 > 2 so that /3 C (^) < 2. For 
the first moment one has ^mi(t) = —f3 + ^.E[L t (0)]. A direct estimate on the local time term 
does not seem useful, but we may examine this equation at discrete times t = 1, 2, ... to find 

mi(n + 1) - mi(n) = -/3 + £7[L„(0)] - £[L„_i(0)]. 

The idea is to show that E[Lt(0)] — E[Lt-i(0)] is bounded below for all t > by a non-zero 
constant. Then the first moment equation shows at least linear growth of the first moment if 
j3 is small enough and hence that f3 c (d) > 0. In the interval [t — l,t] there is probability e~ 2 of 
the event Q, that there is a single jump time r G [i — 1, i — \\ of mass away from the origin, no 
further jumps away from the origin during [t — |,i] and and that no neighbouring sites jump 
onto the origin during [t — 1,1]. Conditional on Q, the local time of the mass at the origin 
evolves, during [r,t], like the local time L for a (rate 2) reflected Brownian motion X with drift 
—(3 starting at zero, that is 

X S = W S - /3s + L s , s > 0. 

Then E[L S ] > E[X S ]. Moroever a comparison argument shows that -E[X S ] is decreasing in f3, 
and when /3 = the variable X s is just modulus of a normal N(0, s) variable. This shows there 
exists ci > 0, independent of /3, so that E[L S ] > c\ for all s G [5, 1], establishing the desired 
lower bound. 



3.1.2 Proof of the asymptotics (4). 

We return to the discrete mass model. Consider a single isolated site undergoing only evap- 
oration at rate p and monomer deposition at rate q = p — /3p 1//2 . Then the rescaled mass 
X t = p~ l / 2 M t is a Markov chain satisfying, on {M t > 1}, 

E[X t+A - X t \a(X s : s < t)\ = f3A + 0(A 2 ), 
E[\X t+A -X t \ 2 \a(X s :s<t)] = (2 - pp~ l ' 2 )A + 0(A 2 ). 

Using these it is straightforward to establish a diffusion approximation: the rescaled mass 
X t = p~ l l 2 M t converge in distribution, as p — > 00, to the diffusion in (16). Extending this to the 
Z d case one expects that, when q = p—j3p 1 / 2 , the rescaled masses {{p- l l 2 M t {x) : x € Z d ) : t > 0) 
converge in distribution to the continuous mass model described above. This convergence 
together with the existence of the critical value /3 c (d), suggest the asymptotics in (1). The 
convergence could be set up, for example, in the space of continuous paths taking values in some 
set of measures on R (simply the empirical measure of the particle positions and masses) with 
a (suitable) weak topology. However convergence for all finite time intervals [0, T] would not be 
sufficient to read off the asymptotic. One would like the convergence of the associated stationary 
distributions (and note the convergence to equilibrium is not expected to be exponential). We 
leave these ideas for later, but we quickly show that the argument for the continuum case applies 
approximately to the (suitably scaled) lattice case, although this inevitably leads to the cruder 
upper and lower constant asymptotics that do not reflect the continuum critical value. Define 
scaled moments by fhk(t) = p~ k ^ 2 mk(t). Choose q = p — ftp 1 / 2 . In these scaled variables the 
the second moment inequality (12) becomes 

fh\{t) - Pmi (t) + (1 - p" 1/2 /3) > for all t > 0. 
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Then if /3 > 2 and p is large enough we see that mi(t), and hence also mi(t), remains bounded. 
This shows that limsup^^ p _1 / 2 (p — q c (j>)) < 2. The scaled first moment equation (9) becomes 

^ = -P +p V*(l- s{t )). (17) 

The argument for the continuum model suggests bounding 1 — s(t) = P[M t (0) = 0] from below. 
For t>l, there is a probability e~ 2 of the event fi that there is single jump time r 6 [t — l,t—\] 
where the mass at the origin jumps to a neighbour, no further jumping of mass from the origin 
during [t— \, t], and no mass jumps into the origin during [t — 1,1]. Conditional on 17, the mass 
at the origin during the interval [r, t] follows a reflected random walk, increasing by 
and decreasing by one at rate p, and starting at zero. For a bound from below of P[Mt(0) = 1 0] 
we may compare, using a simple coupling, with a reflected symmetric random walk M at rate 
p (that is we may, as in the continuous mass case, take (3 = 0). But then exact formulae for 
simple symmetric random walks shows that there exists c x > so that P[M S = 0] > c-jp -1 / 2 for 
all p > 1 and s G [\, 1]. Using this estimate in (17) shows that when /3 is small enough the first 
moment grows at a strictly positive rate, indeed that p _1//2 (<7 c (p) — p) > c^e -2 , for all p > 1, 
completing the crude asymptotics. 



3.2 Small p, q in d = 1 

In d = 1, for small p, q the system is well approximated, on large space-time scales, by the anal- 
ogous continuous space system of massive coalescing Brownian motions on R with evaporation 
and Poisson immigration. We study this model first. 



3.2.1 Continuous space model 

Standard independent Brownian motions on R instantly coalesce upon meeting. Particles have 
masses with values in N, and masses add at coalescence. Also there is deposition of monomers 
at a Poisson rate q, and mass evaporation of each particle, reducing the mass by one, at rate 
p. Let (XI, Ml : i E It) list the positions and masses at time t. We use Qp'^ to denote 
the distribution for this model considered with zero initial conditions. Under diffusive scaling, 
namely X\ = c~ 1 X^ 2f and M\ = M* 2t for c > 0, we obtain a new system of coalescing particles 

Ft N 

with distribution Q c 2 c 3 ■ We will argue below, by similar arguments to the lattice case, that 
the one parameter system, under Qi'^i nas a single critical a c G (0, 00). Then the scaling 
shows the continuum system under Qp,^ has the critical value of qc' N {p) = a c p 3 ^ 2 . Consider 
a system with distribution Q 1 ' a . Since our argument for a critical value of a will only be used 
as a guide for a corresponding lattice argument, we do not give a careful construction of the 
model or proof of the required monotonicity and dependence properties. We may define the 
occupation density and moments by fixing a test function <f> : R — > [0, 00) with J 4> = 1 and 
setting 

s{t) = Ej2HX l t ), m k (t) = Ej2\Mi\ k cj>(Xi). 

i i 

Note by translation invariance that any such test function <f> defines the same quantity. Mono- 
tonicity still implies that s(t) and mkit) (for k > 0) are non-decreasing in t and a and we 
set 

a c = inf{a : mi (00) = 00}. 

The first moment equation for m\(t) is 

dmi , . , 

— (t) = a-s(t). 
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We may bound s(t) from above by s(t) the occupation density for the system with zero evapo- 

Ft N 

ration, that is under the distribution Q Q ' . For this simpler system the scaling, as above, shows 
that the stationary state has a scaling property. Indeed the stationary occupation distribution 
under Qq^ is a scaled copy, by a factor of a 1 / 3 , of the stationary distribution under Qq'i ■ 
Thus s(oo) = C2CK 1 / 3 for some C2 G (0, 00) (in fact C2 is known - see the appendix). Substituting 
in s(oo) < C2CK 1 / 3 into the first moment equation we find that a c < c^J 2 . To find a lower bound 
on a c we consider the second moment 777,2 (i). Again it seems difficult to use a differential in- 
equality (a problem familiar from the density decay estimates for coalescing Brownian motions 
in Bramson and Griffeath [3]). Instead we watch 7772 (t) along integer times t = 1,2,... and 
bound 777,2(77 + 1) — 7772(77) from above. The argument is a little fiddly as each of the dynamics 
(coalescence, evaporation and deposition) contribute, so we start with a heuristic overview. The 
contribution to 7772(77 + 1) — 777,2(77) due to coalescence will be bounded by Cmi(n) 2 by using 
negative correlation for the mass distribution at time 77. The contribution due to the rate one 
evaporation will be bounded by —27771(77) + s(oo) < —27771(77) + C2« 1 ^ 3 . Finally the deposition 
of monomers will contribute 2Ca777i(77) + Co? . Monotonicity of 77 — > 7772(77) leads to 

< 7772(77 + 1) - 7772(77) < Ct77i(t7) 2 + 2Cami(n) - 27771(77) + Co? + c 2 a 1/3 . (18) 

For small a the quadratic inequality prevents 7771(77) from entering the interval (^, ^). But 
the first moment equation shows 7771(77 + 1) —7771(77) < a. This prevents 7771(77) jumping over this 
forbidden interval and forces it to remain bounded for such small a. Hence the critical value 
a c is strictly positive. Now we expand the above argument with a more careful construction 
of the process between times 77 and 77 + 1. The construction exploits the fact that systems of 
coalescing particles can be constructed particle by particle, with each new particle following 
an independent motion until it first hits a previously constructed path, and then following the 
path of the particle it hits. List the (random) positions and masses of the particles at time 77 
as (zi, rrii : i G N). We will give a pathwise construction of the system at time 77 + 1 in three 
steps. 

Step 1: Run CBMs. Run coalescing Brownian paths (independent from the evolution up to 
time 77) starting from (z{ : i G N) over the time interval [0, 1]. 

Step 2: Add deposition. Construct (independently) the rate a Poisson points (yi,U : i G N) 
in R x [0, 1] of the monomer depositions. Then construct (independently) coalescing Brownian 
paths started from these points that also coalesce with any paths in step 1 . Now attach masses 
to the particles by tracking the mass through each coalescence. This yields particle positions 
and masses (Xj, M t * : i G I t ) for t G [0, 1]. 

Step 3: Add evaporation. Add (independently) Poisson marks onto the paths of these co- 
alescing Brownian motions to indicate the times of rate one evaporations. Define the adjusted 
masses (M t l : i G It) along the paths due to these evaporations, (so that some Ml may take 
value 0). 

After all three steps the collection (XI, Ml : i G It, Ml 7^ 0) is then a realization of the true 
system at time n + t. We now estimate the changes in the second moment through each of 
these steps in the construction. Consider the change in the second moment after steps 1 and 2, 
namely 

A 1 :=^0(X{)\Mi\ 2 - X>(^K 2 - 

ieh i 

We rewrite Ai in a different way. For each particle (zi,mi) we let Z % denote the final position 
of this mass at time t = 1 after step 1. Thus Z l will agree with one of the positions (X{ : % G I\) 
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but due to coalescence we may have Z % = Z 3 for many i ^ j. Similarly for each position (yi,U) 
we let Y % denote final position of mass one monomer that was deposed at (yi, U). Then 

Y,^X\)\M[\ 2 = I ^m jX (^=Xi) + ^ X (y fc =Xi) 

i£l ieh V j k 

+ Y J <P(Y k ) X {Y k = Y k ') + 2j2m j( f>(Z j )x(Z j = Y k ) 
k,k' j,k 

=: Ai,i + Ai, 2 + Ai, 3 + A M . 

Each position Zj is, conditionally on <j((xj,mj : i G N)), a N(zj,l) variable. Hence, writing 
(Pt) for the Brownian semigroup, 

P[A 1;1 ] = Ej^PMzjWj = ™ 2 (n) 

by translation invariance. Each pair (Zj,Zk), conditionally on a((xi,rrii : i G N)) and when 
j k, are the positions of a coalescing pair of Brownian motions at time t = 1 started at 
(zj,Zk). We choose 4>(x) = |x(M < 1) an d use the following crude estimate. Suppose (B x ,B y ) 
denotes a pair of coalescing Brownian motions started at (x,si)and (y, s 2 )- Then there exists 
C3 so that for all < si < s 2 < 1 and x, y G R 

P [|Pf I < 1 and B x , B^ coalesce before time 1] < c 3 (|x|" 2 A l)(|x - y\~ 2 A 1). (19) 

See the appendix for a short derivation. Using this we may bound the expectation 

P[A li2 ] < cjfi^mjmHI^r 2 Al)(|zj-^|- 2 A1) 

< c 3 m 2 (n) / (|z|~ 2 Al)(|z-^ / | _2 Al)cZzd2; / 
= 16c 3 m 2 (n). 

Here we have used negative correlation for the second inequality. Next, using (19) for each pair 
of deposed particles (arriving at Poisson a rate), 



A 1) = 16c 3 a 2 , 



and 



£[Ai, 3 ] < c 3 a 2 Jdy J dy 1 J ds x J ds 2 (\y\- 2 A l)(\y - y'\ 
E[A 1A ] < 2c 3 q Jdy J ds J E^mj{\z\~ 2 A l)(|z - y\~ 2 A 1) = 32c 3 am 1 (n). 



j 

Next we look at the change due to evaporation defined by 

a 2 = J>(xj)|M*| 2 -J24>(X{)\M{\ 

ieh ieh 
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t-dt 



The decrease in mass due to evaporation occurs at a Poisson rate along the paths. When a 
mass m at x suffers an evaporation it produces a jump of size <p(x)(l — 2m). Hence 

E[A 2 ] = E f'j2 <P(X$)X(MI > 0)dt - 2E C £ <f>{X\)M, 
h tell Jo «=i? 

= / s{t)dt-2 I mi(n + t)dt 
Jo Jo 

< c 2 a 1/s - 2m 1 (n). 

Combining the estimates on Ai and A2 justifies the difference inequality (18) given at the start 
with the choice C = I6C3. 

3.2.2 Proof of the asymptotics (5). 

We return to the lattice model and choose q = ap 3 ^ 2 for some a > 0. List at time t the positions 
and corresponding masses of all particles as (X\, Ml : i £ I t ). Then the rescaled process 

should converges in distribution as p I to a system of coalescing massive Brownian motions 
with distribution Q 1 ' a as described above. This convergence, together with the critical value 
a c , suggest the asymptotics in (2). However, here we simply aim for crude asymptotics by 
mimicking the continuous argument for the lattice model. Guided by the continuous space 
argument we may bound the occupancy s(t) by the stationary occupancy probability s(oo) for 
the zero evaporation model. For d = 1 models without evaporation a (time reversal) Markov 
duality tool is applicable and one finds s(oo) = 1 — h q (l) where 



h q (x) = E exp (-q X*ds^J , with r = inf{t : Xf = 0}, 



(20) 



and (Xf : t > 0) is a rate 2 continuous time simple random walk on N started at x. Details of 
the duality argument leading to (20) are in the appendix 4.3. Note that h q G [0, 1] solves 

Ah q (x) = qxh q (x) for x G N + = {1, 2, . . .}, 

where A is the discrete Laplacian on N, with boundary conditions h q (0) = 1 and h q (x) — > as 
x — > 00. To see the scaling of h q (l) in q we define A q (x) = h q (xq^ 1 ^ 3 ), so that A q solves 

AA q (x) = xA q (x) for x G <? 1/3 N+, 

where A is the discrete Laplacian on g 1 / 3 N, with boundary conditions A q (0) = 1 and A q (x) — > 
as x —7- 00. Then one can check that sup x |-A 9 (x) — ^4o(x)| — > as q — > 0, where Aq : [0,oo) — > 
[0, 1] is the corresponding C 2 ([0, 00)) Airy function solving A Q (x) = xAq(x). Moreover one can 
check (details are in the appendix 4.4) that the discrete derivative 

g" 1 / 3 (1 - Ml)) = <T 1/3 (A q (0) - A(9 1/3 )) "»> -^o(O) > 0, (21) 

Substituting s(t) < s(oo) = 1 — h q (l) into the first moment equation dm\/dt = q — ps(t) and 
using the above asymptotics leads to 

limsupg c (p)p- 3 / 2 <(-^(0)) 3 / 2 . 
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For the second moment equation we again mimic the argument for the continuous space Brown- 
ian motion case. We consider the second moments m-zit) at the discrete set of times t n = np^ 1 . 
The scaling above suggests that rri2(tp~ 1 ) is 0(p 1 ^ 2 ), and so we aim for a difference inequality 
of the form: there exists C < cc,po > so that for p < po, a < 1 

< m 2 (Wi) - m 2 {t n ) < Cp- 1/2 mi {t n ) 2 + 2Cam 1 {t n ) - 2m 1 (t n ) + C(a 2 + a 1/3 )p 1/2 . (22) 

This implies for small enough a that the values of p~ l / 2 m\{t n ) cannot lie in (^, 2T;)- However 
the first moment equation implies that p _1 / 2 mi(t n+ i) — p _1 / 2 mi(t n ) < a. So for small enough 
a the first moment is bounded. This gives 

liminf q c (p)p~ 3/2 > 0. 
p— »o 

To establish (22) we way repeat the three step construction from section 3.2.1 to construct 
the process between times t n and t n+ \ making the natural changes, namely: replace coalescing 
Brownian motions on R by rate one coalescing random walks on Z; replace rate a deposition 
on [0, 1] x R by rate q deposition on [0, 1] x Z; and replace rate 1 evaporation along the paths 
by rate p evaporation. We may rewite the moments as m^t) = EJ2 X 4'{x)M^{x) for any test 
function <p satisfying 4>{ x ) = 1- We choose <f>(x) = iV~ 1 x(|a:| < p~ 1 / 2 ) where N p 2p~ x l 2 is 
chosen so that the constraint holds. Decompose m,2(t n+ \) — m2(t n ) into terms Ai, A 2 as before, 
where we work however over a time interval of length p^ 1 . We use an analogue of (19), which is 
derived in the appendix 4.2. Let (X x , X y ) be a pair of coalescing rate one simple random walks 
started at (x, Si)and (y, s 2 ). Then there exists c 3 so that for all < si < S2 < p~ l , x,y € Z 
and p < 1 

\Xp-i\ <p _1/2 and X x , X y coalesce by time p' 1 < c 3 {p' l \x\' 2 M){p' l \x-y\' 2 M). (23) 
Following the steps from section 3.2.1, we find E[A\ t i\ = m^im) and 

£[Ai, 2 ] < c^N~ 1 E^mjmj'{p~ 1 \x\~ 2 M^p'^x - y\~ 2 A 1) 

< c z N; l m 2 (t n ) ^(P'Vr 2 A l)(p~\z - z'l 2 A 1) 

z,z' 

< g4^l(*n)p" 1/2 - 

The second inequality uses the negative correlation in (8), and the third inequality uses N p < 
p~ x l 2 and Y1 Z (P~ 1 \ Z \~ 2 ^ -*-) — C ^P ~ 1 ^ 2 f° r an P — 1- I n a similar way we find 

1 / 2 / 

E[A 1>3 ] < -c' 3 cla 2 p 1/2 and E[A 1A ] < -c 3 c\ami{t m ) . 

Using (21) we may choose po < 1 small enough that 

s(oo) < (-2^(0))g 1 / 3 = (-2^(0))a 1 /V /2 - 

Then E[A 2 ] < s(oo) - 2mi(t n ) < (-2A' (0))a 1/3 p 1/2 - 2mi(t n ). Together these estimates imply 
(22) with C = (-2A' (0)) V \c 3 c\, completing the proof. 
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3.3 Remarks on small p, q in d > 3 



There is no continuous space model approximation in d > 1, and the intuition is somewhat 
different. We give a purely heuristic argument suggesting that q c (p) ~ p 2 as p, g 4- 0. The 
argument exploits the small occupation densities and the transience of random walks in d > 3. 
It follows the intuitive argument given in [2] for the modified rate equations that gives the 
decay of the occupation density in coalescing random walks. In [2] they exploit large times to 
obtain low densities, where we exploit small p, q values. We use the informal notation z — > z' 
to denote the event that the particle at z moves to position z' . Choose 0<<r<<iso that 
p-2/d r << p-i choose q = 0(p 2 ). The first moment equation shows that the occupation 
density satisfies s(t) < q/p = 0(p). Then the following approximations seem reasonable for 
small p: 



P[M t (0) > 0, M t (e) > 0] 



P[U a#y {M t _ r (x) > 0, Mt_ r (y) > 0, x -»• 0,y 
Y,PWt-r{x) > 0, M t _ r (y) > 0, x -> 0,y -> 



0] 



^p d p r (x,0)p r (y,e)P[M t _ r (x) > 0, M t _ r (y) > 0] 



^p d p r (x,0)p r (y,0)P[M t _ r (x) > 0]P[M t _ r (y) > 0] 



(24) 



= PdS 2 {t-r) 
~ PdS 2 (t). 

We expect all expressions in these approximations to be 0(p 2 ). The first approximation throws 
away the event that the particles at 0, e were deposited between t — r and t. The expected mass 
of deposited particles during [t — r,t] that reach at time t is of order qr « p, suggesting 
this approximation is valid. The second approximation uses low densities: the error should 
involve three occupied sites and be of order 0(p 3 ). The third approximations involves the 
probability P[x — > 0, y — > e without coalescing]. Evaporation at rate p can be ignored for these 
two particles over this interval since r « p" 1 . Considering the time reversed paths we need to 
find P[0 — > x, e — > y without coalescing], which by transience of the walks is well approximated, 
for \x — y\ » 1, by pdP r (0, x)p r (0, y). Negative correlation shows that the fourth approximation 
is an upper bound. Moreover 

J2Pr(x,0)p r (y,e)P[M t _ r (x) > 0, M t _ r (y) > 0] 



= E 



> E 



(J2Pr(x,0)x(M t - r (x) >0)) 2 

X 

^ Pr (x,0)x(M t _ r (x) >0) 



J2Pr(x,0)P[M t - r (x) >0] 

x 

- s(t - r)J2Pr(x,U) 



= s\t-r)-s(t-r)Y,Pr(x,0). 



The error in the fourth approximation is thus bounded by 

s(t - r) Y,Pr(x, 0) < Cpr~ d / 2 

X 

which is o(p 2 ) since p~ 2 l d « r. Similar approximations are used in [2], where careful error 
bounds are established. We show that just the approximation (24) suggests the scaling for q c (p)- 
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Note 



0<-^-(t) = q(l - s) - pP[M t (0) = 1} - P[M t {0) > 0, M t (e) > 0] 
< q-P[M t (0) >0, M t (e) >0]. 

Using the approximation (24) above suggests that s(oo) < p d 1 ^ 2 q 1 ^ 2 . Using this in the first mo- 
ment equation then suggests that limsupp^g 1c{p)p' 2 < P^P 2 ■ This compliments the rigourous 
upper bound, but we believe we can guess the correct asymptotic. For the case of purely coalesc- 
ing random walks in d > 3, with initially one monomer per site, the occupation density drops 
as p~^ l t~ l . Arratia [1] showed that one can rescale space by t~ l / d to leave an asymptotically 
constant density of occupied sites, and the positions of the rescaled sites converge to a Poisson 
point process on H d . Moroever the mass of the rescaled particles become independent, leading 
to a compound Poisson limit for the empirical measure. The aim is to mimic this for a rescaled 
version of the stationary distribution of our model as p — > 0. There is similar underlying intu- 
ition for a Poisson limit, as follows. Run the process in its stationary distribution for an interval 
of length t. If t « p^ 1 then in bounded regions of space one expects negligible amounts of 
deposition or evaporation. Provided t » p~ 2 l d the particles will move far beyond the typical 
inter-particle distance 0(p 1 ^ d ). Due to transience the process is therefore well approximated 
by free independent motion over this interval, and free motion converges to a Poisson limit. 
To find the approximating compound Poisson approximation we use a more general version of 
approximations such as (24). Take q = ap 2 for some a > 0. We expect 

S[F(M oo (0))F(M oo (x))] « p d {x)E[F{M oo {0i))] E[F(M 00 (x))] (25) 

for all bounded F : N — > R satisfying F(0) = 0, where Pd{x) is the probability of a simple 
random walker never hitting the origin starting at x. Taking F(m) = z m — 1 allow us to 
characterize the mass via its transform, and already suggests that the set of massed particles 
will also converge to a compound Poisson limit. To find the Poisson intensity s a and the 
generating function (j) a {z) for the attached masses, we start by developing ^Ez Mt ^ and taking 
t — > oo to find in the stationarity 

= q(z- l)Ez M ^ + p{z- 1 - l)Ez M ^ X (M t (0) > 0) + 1 - 2Ez M ^ + Ez M °°W +M °°W. 

Define <j>{z) := E^^^M^O) > 0] so that 

Ez M x (o) = Ee- eM °°W X (M t (0) > 0) + P[Moo(0) = 0] = s(oo)0(0) + 1 - s(oo). 

Apply (25) to approximate for small p 

1 - 2Ez M ^ + Ez M °°W +M °°& w p d (l - Ez M ^Y . 
We find, for small p, 

0^q(z- l)(s(oo)0(z) + 1 - s(oo)) + p(z^ - l)s(oo)<j>(z) + Pd s 2 (oo)(l - 4>{z)) 2 . 
Now we suppose, as p — > with q = ap 2 , that 

jO _1 s(oo) ->• S a , <t>(z) ->• 4> a (z) 

and obtain 

= a(z-l) + (z- 1 - l)s a 4> a (z) +p d s 2 a {l - 4> a {z)) 2 . 
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Solving, with the correct sign, we find 



2 Pd s a z(Mz) - 1) = -(1 -z) + v / (l-z)Q(z) 

where Q(z) = 1 — (Ap^asa + l)z + Ap^az 2 . We need to discover the unknown s a . We now 
follow the steps of the mean field analysis in [12] (although our case is slightly simpler in that 
Q is quadratic, due to reducing to a one parameter model). The tail of the variable M with 
generator 4> a can be found by inverting the generating function via the Cauchy integral 

P[M = m} = (2™)- 1 J c ^dz. 

where C is a contour around the origin. The contour can be deformed so that it goes around the 
singularities of <j) a . We expect a growth phase where the first moment predicted by <j) a is infinite 
and s a = a. Solving for the roots of Q(z) when s a = a gives roots z = 1 and z = (Ap^a) -1 . Let 
a c = (Ap^)^ 1 - Then, for a > a c , 4> a has only a single singularity due to the branch cut for the 
root z = (Apda) -1 > 1, and the inversion then shows the tail P[M = m] drops exponentially. 
For a < a c we consider values of s a G [0, a]. Looking at the sum of the roots of Q(z) one finds 
that there must be a root with a real part less than 1 and the inversion would then predict 
exponential growth of P[M = m], except if there is a double root, which would not lead to a 
singularity for cfi a . The unique value for s a leading to a double root is (A^Jp^a — l)/4p d . We 
summarize below our beliefs. 

Conjecture. Under the stationary distribution with q = ap 2 , the scaled empirical measures 
J2 X Moo(x)S p i/d x converge in distribution as p 1 0, in the space of locally finite point measures 
with the vague topology, to a compound Poisson limit ^ i Mi5 Xi , where where (Xi) are Poisson 
with intensity s a and the masses (Mi) are, conditionally, independent identically distributed 
with Laplace transform (f) a (0). 

4 Appendix: details for some proofs. 
4.1 Proof of Proposition 3 and Lemma 4. 

We first consider the comparison theorem and uniqueness for solutions to (6). Suppose M,M 
are two solutions, with the same driving Poisson processes, satisfying Mq < Mq a.s. Then 

d X {M t (x) > M t (x)) 

= -x(M t -(x) = 1, M t -(x) = 0)dP[ p) - X (M t -(x) > M t -(x))J2dPt(x,y) 

+ Y J (x(M t -(x) + M t _(y) > M t -{x) + M t -{v)) ~ x(M t -(x) > M t _(x))) dP t (y,x) 
< J]xW-(»)-^-(») > 0)dPt(y,x) 

y~x 

where in the final inequality we have used x( x + y > 0) — x( x > 0) < x(y > 0)- Then, for 9 > 0, 
Ej2e~ elxl X (M t (x) > M t {x)) < e 9 ( E^2e e ^x(M s (x) > M s {x))ds 

x ^° X 

and Gronwall's inequality implies that Mt < Mi a.s. This implies the comparison theorem and 
pathwise uniqueness of solutions. As usual, uniqueness implies the Markov property. Let M be 
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the solution started from the minimal initial condition Mq = and write v s for the distribution 
of M s . The comparison theorem allows one to run a larger coupled solution M started from 
a random initial condition Mq with distribution v s (but independent of the Poisson drivers). 
Then for non-decreasing F one has 

E [F(M t )} < E [F(M t )] = E [F(M t+s )} 

by the Markov property. This establishes the monotonicity in t in Lemma 4. A simpler argument 
establishes the monotonicity in p and q. Let M be a solution to (6) and let M be a solution 
with the drivers dP^ q \x) replaced by dP iq \x) + dP^~ q \x), where (P^- q \x) : x £ Z d ) is a 
further independent family of rate g — q > Poisson processes. Another Gronwall estimate as 
above shows, if Mq = Mq a.s., that M% < Mt and this establishes the monotonicity in q. From 
zero initial conditions one may check, as before, that E[M t (0) — M t (0)] is non-decreasing, and 
so 

< ^-E[M t (0) - M t (0)] = (q - q) - P P[M t (0) = 0,M t (0) > 0]. 



(it 



Writing sta)(t),s^(*) for the occupancy probabilities of Mt,Mt, we find 



s (5)( t ) _ s W(j) = p[M t (0) = 0,M t (0) > 0] < /p. 

Letting i — >■ oo one deduces, as in section 2.1, that s^(oo) < g/p for s > <J c (p)- Monotonicity 
in p is similar, by constructing the evaporation driver as the sum of two independent families 
pip) ^ = pp( x ) -\-pip-v) ( x j Q f ra tes p and p—p > 0. Existence of solutions can be established by 

the usual iteration argument: set = Mq for all t > and define cadlag adapted processes 
by 

dM{ n \x) = -^2M^ 1 (x)dP t (x,y) + ^2M^- 1 (y)dP t (y,x) 



y~x 



dP ( t q \x) - x{M£T L) {x) > 0) dP K t p, {x). 

Note that d(rj,rj') = J2 X e~ e ^x{"n{ x ) / v'( x )) is a complete metric generating the topology on 
S. Using the simple inequality x( x + y ^ x 1 + y 1 ) — x(x / x') < x(y y') one finds 

sup X (Mi")(x) / M^Hx)) < f X {M { ;j l \x) ± M£r 2 \x))(dPp> +Y J dP s {x,y)) 

s <* ^ y~x 

+ f x {Mt~ l \y) ± Mir 2) (y))^^ s (y,x) 
Jo 



r(n-l) 



(p), 



and hence 



supd(M^ n \M^-^) <(2 + e e ) E\d{M^- l \Mi n - 2 ^) 



ds. 



Now the usual argument shows that M^ converges uniformly over compact time intervals 
to a cadlag adapted limit which is the solution to (6) started at Mq. An argument as for 
the comparison theorem shows the solutions are non-negative. For two solutions M, M' with 
deterministic conditions Mq = rj, M = rf, one has similarly 



E 



d(M t ,M' t )] <C(t,e)d( V ,rf) 
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and this implies the Feller property for the semigroup. The moments (7) are established in 
the same manner, that is by developing dM^(x) and establishing a Gronwall estimate. For the 
attractiveness of the stationary distribution it is convenient to construct an extended system 
that allows particles of infinite mass. Let N U {00} be N with 00 added as a discrete site. At 
the slight risk of confusion we continue to write the equations (6) for this extended system, 
but where we assign natural rules for sites with the value 00. We extend the state space to 
S = (N U {oo}) z . We allow addition, subtraction and comparison on N U {00} via 

00 ± n = 00, 00 + 00 = 00, 00 — oo = 0, n < 00. 

(The operation 00 — 00 will only occur when a particle of mass leaves a site due to a random 
walk step.) The metric d(rj, rj') extends to <S and still generates the product topology. The proof 
of existence, comparison and uniqueness of solutions carries over with only the (small) natural 
changes. The system is a true extension, in that solutions to the original system continue to be 
solutions to the extended system, but which never take the value 00. The comparison theorem 
shows there is maximal solution M with initial condition Mq(x) = 00 for all x. Define 

A t = {x e Z d : M t (x) = 00}. 

Examining the evolution of dx(M t (x) = 00) shows that the evolution of A t is precisely that 
of the set of occupied sites for instantly coalescing random walkers, started from all sites in 
Z d occupied. In particular P(x € At) is constant in x and converges to as t — > 00. If M is 
the minimal solution started from Mq = 0, then comparison implies Mt < Mt for all t. We 
claim that Mt{x) = Mt(x) for x £ A\. The proof is another simple Gronwall estimate showing 
E\E xeA f e- \ x \ X (M t (x) < M t (x))\ = 0, by using 

d x (M t (x) < M t (x) < 00) < Y,x(M t ^(y) < M t -{y) < oo)dP t (y,x). 

y^x 

This gives all the ingredients needed to complete the proof of Proposition 3. Starting from zero 
initial conditions, the monotonicity of the Laplace functional implies convergence in distribution 
to Mqo £ S. But the coupling above by M and the decay of P[x 6 A t ] implies that 6 S 
almost surely. The limit distribution is a stationary distribution (an argument that uses 
the Feller property). Any other solution can be coupled between M and M and so converges 
in distribution to v. In particular the stationary distribution is unique. 

4.2 Proof of the coalescent pair estimates (19) and (23). 

We give the argument for the discrete walk case (23), and note the Brownian case can be treated 
entirely similarly. We do not seek accurate bounds. Split the problem in two by using Cauchy 
Schwarz, namely 

(-PO-Xp-il < P~ 1/2 and X x , XV coalesce by time p- 1 ]^ 

< P[\Xp-i I < jT 1/2 ] P [X x , Xy coalesce by time p- 1 } . (26) 

For \x\ > 2p~ 1 / 2 and p<lwe have 

P[\X x p -,\ < p- 1 ' 2 } < P[|X p °-i| > \x\/2] < 16\x\- 4 E[\X^\ = m\x\-\p~ l + Zp- 2 ). 

Thus, for suitable C, PWX*^ < p- 1 / 2 ) < C(p~ 2 \x\- 4 A 1) for all x and p < 1. Take walks 
X z and X y starting at (z, S2) and (y, S2). Write Q x for the law of a simple random walk on Z 
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started at x on pathspace. 



P[X Z , X y coalesce by time p^ 1 ] 

= Qz-y[X hits zero before 2(p _1 — s 2 )\ 

< Q\ z - y \ [X hits zero before 2p _1 ] 

< 2Q\ z - y \ [X 2p -i < 0] by the reflection principle 

< 2mm{Q z _ y [X 2p - 1 < 0],Q y _ z [X 2p -i < 0]}. 

By conditioning at time s 2 we have 

P[X X , X v coalesce by time p^ 1 ] 

< 2E mm{Q xs ^ y [X 2p -i < 0],Q y -x Sa [*2 P -i < 0]} 

< 2 umi{EQ x ^y[X 2p -, < 0],EQ y -x fa [X 2p -i < 0]} 

= ^Q\x~y\[ X 2p- 1 +s 2 ~s 1 < 0] 

< C{\x-y\~ A p- 2 M) 

by a Markov inequality as before. Substituting these estimates into (26) completes the proof. 
4.3 Proof of the duality (20). 

The Brownian web and the dual Brownian web (see [8]) is a useful construction for the study 
of the model in d = 1, at least with zero evaporation. Informally, a realization of the web 
on [0,i] x R has coalescing Brownian paths starting from every space-time point (t,x). If one 
superimposes an independent Poisson a process on [0, t] x R, yielding the space-time points 
(ti, Xj)j 6 N of the deposed monomers, then almost surely the web has a unique set of coalescing 
paths starting at each (ti,Xi). These can be used to construct the process over [Q,t] under 
the distribution Q Q ' a . The dual web contains an additional set of coalescing paths running 
backwards over the interval [Q,t], starting at every point in {t} x Q. Moreover, almost surely, 
a forwards path, started from (ti,Xi), never crosses any of the backwards dual paths. This 
property implies the following identity for the probability, under Q ' a , of an empty interval 
[x, y] at time t, where x < y, is given by 

P [[x, y] is empty at time t] = E [[B% , B y s ] n (U,Xi) ieN = for all s G [0, t]\ 

where Bg , B y s are the positions at time s of the coalescing backwards dual paths started at 
Bf = x and B\ = y. The point is that if any particle is deposited in the interval (B*, B y ) then 
its forwards path would have to end up in the interval [x, y] at time t. Taking expectations over 
the Poisson deposition yields 



P[[x,y] is empty at time t] = Eex.p 




Letting t — > 00 one finds the empty interval probability for the stationary distribution under 

Qo, a 1S s iven b y 

P [[x,y] is empty at time 00] = E exp ^—a J X y ~ x ds^J . 

where X v ~ x is a rate 2 Brownian motion started at y — x and r = inf{s : X y ~ x = 0}. The right 
hand side is given by the Airy function Ao(a 1 ^ 3 (y — x)). Differentiating in x and then letting 



23 



y 1 x gives the stationary occupation density as — q 1 / 2 ^4q(0). This identifies the constant C2 
in section 3.1 as C2 = -Aq(0). One can see many related arguments in [21]. A more detailed 
description of the evolution and stationary distribution under Q Q ' a is in preparation in a current 
thesis. A continuous time random walk analogue of some of the above construction is possible, 
quite close to the graphical construction of the voter model and its dual. Start with a graphical 
construction of coalescing rate one simple random walks on the lattice [0, t] x Z as follows. Let 
(U, yi, ±)j 6 N be the points of a rate one marked Poisson process on the dual lattice [0, t] x (Z+^). 
The marks ± are chosen independently for each point, with equal probability. Passing a point 
(t, n+ +) a forwards path at n will jump to n + 1. Passing a point (t, n+ g, — ) a forwards path 
at n + 1 will jump to n. These rules define a unique right continuous forwards paths started at 
any point in [0, t] x Z. Moreover paths started at different points will coalesce instantaneously 
should they meet. The same marked points can be used to define backwards paths living on the 
dual lattice. When such a path approaches a point (t, n+ |, +) it will jump to n — |, and when 
it approaches a point (t, n + 3, — ) it will jump to n + \. These rules define right continuous 

paths (Y s 2 : x G Z, s G [0,t]) so that their backwards versions s — > Y t x _ s are coalescing rate 

x+ i , 

one random walks started at X t 2 = x + ^. One can then check that forwards paths never 
jump over backwards paths in the following sense: if a forwards path X and a backwards path 
Y are both defined at s and t, then 



X s < Y s if and only if X t <Y t . 



Now the argument is as for the Brownian web. Lay down an independent set of deposition 
points (ti, Xj)i £ N as a rate q Poisson process on [0, t] x Z, and use the forwards paths constructed 
above to build the process on [0, t] that has distribution Q ' .A site x G Z is empty at time t 

precisely if the interval [Y s 2 ,Y S 2 ] contains none of the deposition points (ii,Xi)ieN at any 
time s G [0,i\. Now mimic the Brownian argument to reach (20). 



4.4 Proof of convergence to the Airy function and (21). 

Note that A q and Aq are convex, non-increasing and take values in [0,1]. Let DA q (x) = 
q~ 1 / 3 (A q (x+q 1 ^ 3 ) — A q (x)) be the discrete derivative, so that the discrete Laplacian is AA q (x) = 
q-V^DAgW-DAqix-q 1 / 3 )) onq 1 / 3 N + . Using A q (y) -A q (x) = q 1/3 Y, x < z<y DA q (z), and the 
fact that x — > DA q (x) is non-increasing, one finds that \DA q (x)\ < x~ l for all x and q. Together 
with the identity /S.A q {x) = xA q (x), this shows that both DA q (x) and AA q are bounded over 
intervals [0, K] uniformly in q. This gives the required equicontinuity to deduce that there is a 
continuously differentiable A : [0, 00) — > [0, 1] so that, at least subsequentially, 



max \A a (x) — A(x) 

x<K q 



and max|ZM g (x) - A'{x)\ -)• 0. 

x<K 



We aim to check the limit points A are unique and solve the Airy equation. Passing to the limit 



m 



tV3 V 



x<z<.y 



zA q (z) 



q 1 ' 3 V AA q (z) 



x<z<y 



leads to J* 00 zA(z) dz < x 1 and hence that A{z 
A(0) 



< —DA q (x) < x~ l 
as z — > 00. The other boundary condition 



1 is immediate. Another passage to the limit in 

1 + gVs £ (D q (y) + q^ £ zA q (z)) 

0<y<x 0<z<y 



A q (x) 
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gives A(x) = 1 + Jq(A'(0)+ J Q y xA(z)dz)dy, completing the proof that A 6 C 2 ([0, oo)), a solution 
to the Airy equation. Uniqueness of solutions is straightforward via a maximum principle 
argument. 
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